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ABSTRACT 

A detailed study is presented of the decaying solar active region NOAA 10103 ob- 
served with the Coronal Diagnostic Spectrometer (CDS), the Michelson Doppler Im- 
ager (MDI) and the Extreme-ultraviolet Imaging Telescope (EIT) onboard the So- 
lar and Heliospheric Observatory (SOHO). Electron density maps formed using Si 
x(356.03 A/347.41 A) show that the density varies from ~10^^ cm~^ in the active re- 
gion core, to ~7x 10^ cm~^ at the region boundaries. Over the five days of observations, 
the average electron density fell by ~30 per cent. Temperature maps formed using Fe 
xvi(335.41 A)/Fe xiv(334.18 A) show electron temperatures of -2.34x10^ K in the 
active region core, and ~2.10x 10^ K at the region boundaries. Similarly to the electron 
density, there was a small decrease in the average electron temperature over the five 
day period. The radiative, conductive, and mass fiow losses were calculated and used 
to determine the resultant heating rate {Ph)- Radiative losses were found to dominate 
the active region cooling process. As the region decayed, the heating rate decreased by 
almost a factor of five between the first and last day of observations. The heating rate 
was then compared to the total unsigned magnetic fiux {^tot = / d^\Bz\)^ yielding a 
power-law of the form Ph ~ ^^'^^^^^'^'^ . This result suggests that waves rather than 
nanofiares may be the dominant heating mechanism in this active region. 
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1 INTRODUCTION 

Since the discovery of highly ionized species of iron in the so- 
lar corona in the 1930's, physicists have been puzzled by the 
high temper atures observed in the outer solar atmosphere 
(lEdlenlll937l) . It is widely accepted that the magnetic field 
plays a fundamental role in the heating process, but precise 
measurements of the coronal magnetic field are currently 
impossible. Indirect methods are therefore adopted which 
rely on measurable quantities such as the electron temper- 
ature, density and photospheric magnetic flux. These mea- 
surements can then be compared to theoretical predictions. 

Models for coronal heating typically belong to one of 
two broad categories. In wave (AC) heating, the large-scale 
magnetic field essentially acts as a conduit for small-scale, 
high-frequency Alfven waves propagating into the corona. 
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For constant Alfven wave amplitude (v^), the total power 
dissipated in an active region is, 

Ph = \f^{v^)^tot ergs s-\ (1) 

where p is the mass density, and ^tot is the total unsigned 
magnetic fiux, 

^tot = J dA\B,\ Mx, (2) 

where Bz is the longitudinal component of the magnetic 
field. 

In stress (DC) heating, the coronal magnetic field stores 
energy in the form of electric curr ents u ntil it can be dissi- 
pated, e.g., by nanofiares (ParkelEoS^ . The total power 
can be estimated by, 

PH--\v\^tot ergss-\ (3) 
and the constant of proportionality describes the eflaciency 
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Figure 1. The evolution of NOAA 10103 during September 10—14. Top row: EIT (195 A) images. Middle row: MDI magnetograms. 
Bottom row: MDI continuum images. The inner boxes in each image denote the field-of-view of the CDS raster. 



of magnetic dissipati on, wh ich might involve the random 
footpoint velocitv. v (|Parkei| |l983,') . or simply the geom etry 
(jBrowning. Sakuria. Pries tlll98 rtlFisher et al■lll998^ . 



Several authors have linked the photospheri c magnetic 
flux to EUV and X-ray line intensity. Gurm an et al.l (11974) 
found that the hne intensity of Mg x (624.9£A) was pro- 
portional to the magnetic flux density. Schrlive? ("1987*) re- 
lated the integrated intensities of chromospheric, transi- 
tion region, and coronal lines to the total magnetic flux 
by a power law, the index of which was depe ndent on the 
scale h eight. This result was later confirmed bv lFisher et al.l 
(|l99 8^. who showed that X-ray luminosity is highly cor- 
related with the total unsigned magnetic flux. Van Driel- 
Gesztelyi et al. (2003) also showed a power-law relationship 
between the mean X-ray flux, temperature, and emission 
measure, and the mean magnetic field by studying the long- 
term evolution of an active region over several rotations, at 
times when there were no significant brightenings. 



In this paper, we study the evolution of a decaying ac- 
tive region using the diagnostic capabilities of the Coronal 
Diagnostic Spectrometer (CDS; Harrison et al. 1995) and 
Michelson Doppler Imager (MDI; Scherrer et al. 1995) on- 
board SOHO. Using the temperatures, densities, and dimen- 
sions of the active region, the heating rate is calculated and 
compared to the total unsigned magnetic fiux. These results 
can be put in the context of theoretical models. Section |21 
gives a brief overview of the active region, a summary of the 
instruments involved, and a description of the data analysis 
techniques. Our results are given in Sectional and discussion 
and conclusions in Section HI 



2 OBSERVATIONS AND DATA ANALYSIS 

NOAA 10103^ was observed by CDS, EIT, and MDI for five 
consecutive days during 2002 September 10-14. Fig.Qshows 
the general evolution of the region over that period. 



2.1 The Coronal Diagnostic Spectrometer 

EUV spectra were obtained with the CDS instrument, which 
is a dual spectrometer that can be used to obtain images 
with a spatial resolution of ^8 arcsec. The Normal Inci- 
dence Spectrometer (NIS), used in this study, is a stigmatic 
spectrometer which forms images by moving the solar im- 
age across the slit using a scan mirror. The spectral ranges 
of NIS (308-381 A and 513-633 A) include emission lines 
formed over a wide range of temperatures, from '^10^ K at 
the upper chromosphere, through the transition region, to 
'-^lO^ K at the corona. The details of the AZ_ddep1 observ- 
ing sequence used in this study can be found in Table Q] 

The raw CDS data were cleaned to remove cosmic rays, 
and calibrated to remove the CCD readout bias and convert 
the data into physical units of photons cm~^ s~^ arcsec"^. 
Due to the broadened nature of post-recovery CDS spectra, 
the emission lin es were fitted with modified Gaussian profiles 
as described bv iThompsonl (|l999[ ). The Gaussian term was 
defined as. 



G{\) — exp 



(4) 



and the wings by. 



1 See http: //www.solarmonitor.org/20020910/0103.html| 
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Figure 2. Portions of the CDS spectra from NOAA 10103 on September 10, together with backgrounds (horizontal dot-dash hnes) and 
broadened Gaussian fits. 



Table 1. Details of the CDS AZ_ddep1 observing sequence 



Parameter 



Value 



Date 

Region name 
Instrument 

Wavelength range (A) 
Slit size (arcsec-^) 
Area imaged (arcsec^) 
Exposure time (s) 
Number of slit positions 



2002 September 10-14 
NOAA 10103 
CDS/NISI 

332-368 
4.064x240 
243.84x240 
50 
60 



Table 2. Ions, Formation Temperatures, Transitions and Wave- 
lengths of Emission Lines Identified in This Work 



Ion 


Log Te 


Transition 




A/A 


Si X 


6.1 


2s22p2 3p^/2-2s2p 


^3/2 


347.41 


Si X 


6.1 


2s22p 3p3/2-2s2p2 


^1^3/2, 5/2 


356.03 


Fe XIV 


6.3 


3s23p 2p^/2-3s3p2 


^D3/2 


334.18 


Fe XVI 


6.4 


3s 2Si/2-3p 2p3/2 


335.41 


Fe XVII 


6.7 


2p^ 3s 3pi-2p5 3p 


'B2 


347.85 



W{X) = 



(^) +1' 



(5) 



where A is the wavelength, Ao is the central wavelength of 
the line, a is the Gaussian width, and is the FWHM. 

The combined function describing the line profile can 
be expressed as. 



B{X) ^ - a)G{X) + aW{X)], 



(6) 



where / is the amplitude of the line profile and a can be the 
amplitude of the red or the blue wing. The line flux is then 
given by. 



/ 



B{X)dX 



dred 



V Cired / 



27V 



(7) 



These broadened Gaussian profiles were then fitted to 
emission lines in three wavelength intervals using the xcfit 
routine in the CDS branch of the SolarSoftWare tree (SSW; 
Freeland & Handy 1998; see Fig. 12)). Two of the intervals 
were centreed on each of the density sensitive Sl x (347.41 A) 
and Sl X (356.03 A) lines, while the third contained the tem- 
perature sensitive Fe xiv (334.18 A) and Fe xvi (335.41 A) 
pair. Due to the relatively low intensity of Fe xiv (334.18 A) 
compared to that of the adjacent Fe xvi (335.41 A) transi- 
tion, an upper constraint on the width of cr = 0.4 A (FWHM 
= 0.17 A) was placed on the Fe xiv hne. The primary lines, 
their formation temperatures, transitions, and rest wave- 
lengths are given in Table 



2.2 The Michelson Doppler Interferometer 

Magnetic field measurements were taken by the MDI instru- 
ment, which images the Sun on a 1024 x 1024 pixel CCD 
camera through a series of increasingly narrow filters. The 
final elements, a pair of tunable Michelson interferometers, 
enable MDI to record filtergrams with a FWHM bandwidth 
of 94 mA. Several times each day, polarizers are inserted 
to measure the line-of-sight magnetic field. In this paper, 
5- minute- averaged magnetograms of the full disk were used, 
with a 96 minute cadence and a pixel size of 2 arcsec. 

iBerger LiteJ (|2003l) analyzed Advanced Stokes Po- 
larimeter (ASP) and MDI magnetograms, and found that 
MDI underestimates the magnetic flux densities by a factor 
of 1.45 for values below 1200 G. For flux densities higher 
than 1200 G, this underestimation becomes nonlinear, with 
the MDI fluxes saturating at -1300 G. Values below 1200 G 
were therefore corrected by multiplying by 1.45, and values 
above 1200 G wer e approximately co rrected by multiplying 
by a factor of 1.9 JGreen et al.ll2nn3h . 

Before calculating the total unsigned magnetic flux, 
two final corrections were applied to the data. The first re- 
sults from the fact that the measured line-of-sight flux devi- 
ates more and more from a radial measurement as one ap- 
proaches the limb. For simplicity, we assume that magnetic 
fields in the photosphere are predominantly radial. Then, the 
radial field strength becomes equal to the line-of-sight field 
strength times 1/ cos^, where is the heliocentric distance 
of the region from Sun centre. 

Active region areas. A, were calculated by counting all 
pixels above 500 G and multiplying by the appropriate factor 
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Figure 3. The cosine-corrected active region area, total magnetic 
flux density, and total magnetic flux for September 10-14. 



to obtain the active region area in cm^ . This threshold value 
for the magnetic field was found to adequately separate ac- 
tive region structures from neighbouring areas of quiet-Sun 
and plage. As the region approached the limb, the effects 
of foreshortening became significant. Measured areas were 
therefore corrected by dividing by cos^. The resulting total 
unsigned flux was then calculated using Equation 

Modern high resolution space based instruments, such 
as the Transition Region and Coronal Explo rer (TRACE), 
have been used bv lAschwanden et al.l (|200lh to show that 
coronal loops observed in different bandpasses are not nec- 
essarily cospatial. Due to the coarse resolution of CDS, we 
make the assumption that active regions occupy a hemi- 
spherical volume, y = 2/3 tt^^^^A^^^, with a mean loop 
length of L = VttA. 



3 RESULTS 

3.1 Morphology and Magnetic Field 

The top row of Fig. Q shows a series of 360 x360 arcsec 
FIT images obtained in the 195 A bandpass. The first three 
images show a number of loops to the south and north of 
the region, which are not visible from 2002 September 13 
onwards. On September 10, the MDI magnetogram shows 
a simple bipolar p region, which is then classified as a 
on the following day. The region subsequently decreased in 
both size and complexity as it approached the west limb. 
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Figure 4. Theoretical line ratio as a function of density from 
CHIANTI for Si X (356.03 A/347.41 A). The vertical dotted lines 
denote the typical range of densities expected to be found in an 
active region. 



The overall decay of the active region is clearest in the MDI 
continuum images in the bottom row of Fig. □ 

Table^land Fig.^lshow the decay of the region in terms 
of the cosine-corrected area and magnetic flux for Septem- 
ber 10-14. The region was observed to have an initial area 
of 1.26x10^^ cm^ which fell to 6.54x10^^ cm^ by Septem- 
ber 14. The total unsigned magnetic flux density also shows 
a similar trend, falling from close to 4.80x10^ Mx cm^ to 
1.40x10^ Mx cm . The product of the region area times 
the total magnetic flux density is then given in the bottom 
panel of Fig. El As the region decays, the total unsigned 
magnetic flux falls off by a factor of 5-6 over the five days 
from September 10 to 14. 



3.2 Electron Densities 

Flectron density maps were generated using the Si x 
(356.03 A/347.41 A) ratio in conjunction with theoretical 
data from the chianti v4.2 atomic database (Der e et al.l 
Il997f) . Fig. 21 shows a plot of the density sensitive Si X 
(356.03 A/347.40 A) ratio together with the expected range 
of densities in an active region of this size and class. 

The density maps, presented in the second row of Fig. m 
show values of ^10^° cm~^ in the active region core, and 
^^7x10^ cm~^ in the region boundaries. These are in good 
agreement with previous active region density measurements 
(Gallagher et al. 2001; Warren k Winebarger 2003). Densi- 
ties in the core were observed to fall from ^10^^ on Septem- 
ber 10 to 3.9 X 10^ cm~^ on September 14, excluding the high 
density limit ('-^2.5x10^° cm~^) that was reached during a 
C- class flare on September 13. This behaviour is more ev- 
ident in the top panel of Fig. El which clearly shows that 
the average electron density changed by ^30 per cent over 
the five days. Subsequently, plasma that was found to reach 
the high density limit was excluded from any further calcu- 
lations. The average electron densities are listed in Table El 
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Figure 5. Top row: CDS images of NOAA 10103 over the five days of observations. Second row: Electron density maps derived from 
the Si X (356.03 A/347.41 A) ratio. Third row: Electron temperature maps found from the Fe xvi (335.41 A)/Fe xiv (334.19 A) ratio. 
Bottom row: Electron pressure maps from the product of the density maps with the formation temperature of the Si x lines. The panels 
for September 13 show the compact brightening. 



Table 3. Principal Parameters Derived for NOAA 10103 for 2002 September 10-14 



Date 


(Ne) 

lO^cm-^ 


(Te) 

lO^K 


{Pe) 

lO^^K cm-3 


A 

10i8cm2 


V 

1028cm3 


L 

lO^cm 


2002 September 10 
2002 September 11 
2002 September 12 
2002 September 13 
2002 September 14 


2.46±0.31 
2.71±0.53 
3.04±0.79 
2.46±0.52 
1.87±0.20 


2.26±0.03 
2.24±0.03 
2.26±0.03 
2.24±0.02 
2.16±0.02 


3.10±0.60 
3.41±0.84 
3.83±1.16 
3.10±0.65 
2.35±0.43 


12.66±0.75 
10.23±0.64 
10.58±0.76 
9.89±1.03 
6.54±1.22 


1.69±0.15 
1.23±0.11 
1.29±0.14 
1.17±0.18 
0.62±0.17 


6.31±0.37 
5.66±0.35 
5.76±0.41 
5.57±0.58 
4.53±0.84 


Date 


Pr 

lO^^ergs s~^ 


Pc 

lO^^erg s-i 


Pf 

lO^^erg s-i 


Ph 

lO^^ergs s~^ 


^tot 

10^4 Mx 




2002 September 10 
2002 September 11 
2002 September 12 
2002 September 13 
2002 September 14 


3.37±0.90 
2.97±1.12 
3.92±2.07 
2.33±1.05 
0.71±0.24 


0.28±0.07 
0.19±0.05 
0.21±0.06 
0.30±0.07 
0.07±0.03 


0.25±0.04 
0.37±0.08 
0.43±0.12 
0.42±0.14 
0.18±0.06 


3.90±1.57 
3.53±1.92 
4.56±3.10 
2.59±1.40 
0.97±0.68 


5.97±0.35 
4.01±0.25 
3.97±0.28 
2.78±0.29 
0.91±0.17 
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Figure 6. Variations in electron density (A^e), temperature (Te), 
and pressure (Pe) for September 10-14. 



3.3 Electron Temperatures 

E lectron temperature maps were created using the method 
of lBrosius et al.l (|l99(t ) , which employs the ratio of the fluxes 
of various ionization stages of iron. This method obtains 
a polynomial fit to the logarithm of the temperature as a 
function of the logarithm of the emissivity ratio for selected 
line pairs, under the assumption of an isothermal plasma, 
with the form. 



Te = ao + ai(log R) + a2{log Kf + asilog R)^ , 



(8) 



where the parameters, ap, ai, a2 , and as were initially tab- 
ulated bv lBrosius et al.l ([1996^, and R is the intensity ratio 
Fe XVI (335.41 A)/Fe xiv (334.18 A). 

Using the most recent theoretical atomic data from CHI- 
ANTI v4.2, the values of ao,.--? ^^3 were found to change some- 
what and resulted in a slightly higher temperature (^5 per 
cent) than those predicted using'Brosius et all ([1996') values. 
Both previous and updated parameters have been included 
in Table 0] while the resulting temperatures are presented 
in Fig.O In addition, the Fe xiv (334.18 A) line is density 
sensitive, and was accounted for by determining the temper- 
ature across the active region at the corresponding density. 

The temperature maps are displayed in the third row 
of Fig. |51 As expected, the temperature maps show a close 
spatial correlation with the intensity and density maps. As 
the region evolves, the temperature remains constant at 
^2.25x10^ K, with just a slight decrease (by ^4 per cent) 
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Figure 7. The temperature sensiti ve curves for Fe xvi 
(335.41 A)/Fe xiv (334.18 A) for the Br osius et al.l(ll99(j) method 
and the updated CHIANTI v4.2 data, both plotted for plasma den- 
sities of 10^ and 10^° cm"^. 
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Figure 8. Portion of the spectrum from the compact brightening 
on September 13 showing the Fe xvil (347.85 A) emission line 
and neighbouring lines. Fe xvil has a formation temperature of 
- 5 X 10^ K. 



on the final day of observation. More statistically significant 
is the variation of electron temperature across the region, 
which ranges from ^-^2. 10x10^ K in the active region bound- 
ary to ^2.34x10^ K in the core. The values derived for the 
temperature are heavily dependent on the ao coefficient in 
Equation @, which accounts for the restricted temperature 
sensitivity of this method. The temperature is also affected 
by the brightening on September 13 as can be seen from the 
corresponding map. Average temperature values are listed 
in Table El and are plotted in Fig. |S1 Aga in, these results 
are in good agreement with those found in lGallagher et all 
(|200J) using similar methods. 

The C-class flare on September 13 is not only evident 
in the intensity and density maps of September 13, but was 
also identifled due to the presence of the Fe xvii (347.85 A) 
line which has a formation temperature of ~5xl0^ K. A 
portion of the spectrum for the flare is shown in Fig. |S] 
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Table 4. Comparison between the parameters used in Eauationl^ 
for Fe XVI (335.41 A)/Fe xiv (334.18 A). 



Eh — Er + Ec + Eb 



ergs cm ^ s 



(9) 



where Er is the total radiative losses of the plasma, Ec 
is the thermal conductive flux, Ep the energy lost due 
to mass flows, and Eh is the energy required to balance 
these losses (e.g. Antiochos & Sturrock 1982, Bradshaw & 
Mason 2003). The radiative loss term, Er^ in Equation © 
can be written in terms of the electron density, A^e, and the 
radiative loss function, A(Te), 



Er = -NlK{Te) ergs cm"^ s~ 



(10) 



The radiative loss function is usually ap proximated by an - 
alytical expressions of the form N^x'^e (ICook et al.l 
While this method is simple to implement computationally, 
it does not capture accurately the fine-scale structure of the 
radiative loss function. A more appropriate technique is used 
here, which relies on interpolating the valu es obtained from 
CHIANTI, using the coronal abundances of ^Feld man| (|l992|) 
and t he ionization balance calculations of Maz zotta et alJ 
(|l9 98\ The choi ce of coro nal abundances was motivated by 
recent work by Del Zanna Mason. (,2003) who show that 
low FIP elements in active regions have abundances which 
are coronal. 

Assuming classical heat conduction along the magnetic 
field, the conductive flux can be expressed as. 



Ec = 



d_ 

ds 

where the 



i5/2 



dTe 
ds 



ergs cm 



(11) 



thermal 

-1 ^^-1 



con ductivity is = 0.92 x 

10-^ ergs cm-' R-^/^ ( Spit zer 1962"), Te is the electron 
temperature. Various approximatio ns for Equation (IHI) have 
been used in previous studies: fAschwanden et al.l f 1999^ 1 
used an expression that was heavily dependent upon the 
temperature gr adient, (dT/ds)^, in the study of active region 
loops, whereas IVaradv et al.l (|2000l) used an approximation 



Study 


Log Ne 


ao 


ai 


a2 


as 


Brosius et al. C1996) 


9.0 


6.237 


0.111 


0.0087 


0.0016 


10.0 


6.217 


0.108 


0.0076 


0.0018 


Updated values 


9.0 


6.259 


0.111 


0.0105 


0.0015 




10.0 


6.240 


0.108 


0.0097 


0.0015 



Radiative Losses 



026. 



3.4 Electron Pressures 

The pressure maps in the bottom row of Fig. El follow a 
similiar behaviour to that seen in the density. Values for the 
pressure in the high density core of the active region remain 
-^10^^ K cm"^ from September 10-13 and drop to 4x10^^ K 
cm~^ on September 14, again not taking into account the 
effects of the brightening. Similar to the average electron 
density, the average pressure varied by '-^30 per cent between 
September 10-14 (see the values presented in Table and 
Fig. El . 



3.5 Power Balance 

The steady-state energetics of a coronal loop can be ex- 
pressed as. 



~TV2/ds2 

~T5/2 (dT/ds) 
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Figure 9. The radiative, conductive, and mass losses, and the 
heating rate for NOAA 10103 for September 10-14. 



strongly dependent upon temperature, T^^^, for post-flare 
loops. Here, we assume a semi-circular loop geometry and 
approximate Equation (fTT)) by. 



Ec 



5/2 



ATe 

(L/2)2 



ergs cm 



(12) 



where ATe is the difference in the maximum and minimum 
temperatures given in Section lT^ and ds is approximated by 
L/2. This approximation is therefore not overly dependent 
on temperature or the temperature gra dient, but acts as a 
balance be tween the approxini ations of lAschwanden et al.l 
Jl99d) andVaradv et alj J200(tl . 

The energy losses due to mass motions of the plasma, 
Ef, can be expressed as a sum of the kinetic energy and the 
internal energy of the plasma. 



^ ds 



V (^^PV^ + ^NeksTe 



ergs cm 



(13) 



where v is the flow velocity of the plasma, and ds is again 
approximated to be the loop half-length {L/2). Line-of-sight 
velocities at the core of the active region were found to be 
'-^lO km s~^ which is consistant with*Brvnildsen et ajj (^O^) 
who detected flow speeds of ^-^15 km s~^ also using CDS. Us- 
ing these velocity values, as well as typical parameters from 
Table El the mass loss term was found to be ^-^10 ergs s 
Energy losses due to mass motions can therefore be consid- 
ered negligible compared to radiative losses. 
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In orde r to compare the work desc r ibed in this paper 
with that of lFisher et all dTgol) , |Parke3 (1198.-^) . and others, 
Equation Q has been integrated over the volume of the 
active region, and expressed as a power balance equation: 



Ph^Pr + Pc + Pf ergs s" 



(14) 



where Pr, Pc, and Pf are the power lost due to radia- 
tion, conduction, and mass flows, respectively. Ph is then 
the heating rate required to balance these losses. 

Using the active region properties detailed in Table El 
the radiative, conductive, and mass flow losses were then 
calculated using Equations (fT?7|) . (fT^ . and (fT^ . These re- 
sults, together with the heating rate calculated using Equa- 
tion dl, are presented in Tableland Fig. El Due to the 
difficulty in correlating emission seen in CDS with a par- 
ticular magnetic flux concentration, the region's average 
properties were used for comparision with MDI. The sec- 
ond panel down of Fig. El also shows how the conductive 
flux values depend on how Equation is approximated. 
In each case, the conductive losses are found to be much less 
than the radiative losses. Indeed, the average radiative losses 
{(Pr) ~ 3 X 10^^ ergs s""*^) are found to exceed both the con- 
ductive losses {{Pc) ^ 2 X 10^^ ergs s~"^) and the mass flow 
losses {{Pf) ^ 4 x 10^^ ergs s~-^) by approximately an order 
of magnitude. lAschwanden et al.l pOOO) also found an order 
of magnitude difference between the radiative and conduc- 
tive losses in coronal loops despite the actual values being 
somewhat lower than those found here due to the insensi- 
tivity of EIT filter ratio techniques. As can be seen from 
Fig. El the heating rate falls by close to a factor of 5 be- 
tween September 10 and 14 and is not significantly affected 
by whichever conductive loss equation is used. 

The top panel of Fig. E3 shows the heating rate as a 
function of the total unsigned magnetic flux. A least-squares 
fit to the non- flaring data (i.e., neglecting the high density 
plasma of September 13) yielded a power-law of the form 
Ph ^tot^ where 7 = 0.81 ± 0.32. Power-laws with slopes 
of 1 and 2 from Equations (^3 and Q are also shown for 
comparison. The bottom panel of Fig. 1101 shows how this 
relationship varies when densities from the core and the 
boundaries of the active region as well as from the quite 
solar corona (see Fig. El are used instead of average values. 
This wide range of A^e shows that 7 varies from 0.2 in the 
quiet-Sun to 1.5 at the core of the region. It can therefore 
be concluded that average values of these parameters are a 
reasonable representation of the entire active region. 



4 DISCUSSION AND CONCLUSIONS 

A detailed study of an evolving active region has been de- 
scribed using measurements from several instruments on- 
board SOHO. The region was observed to decay in size and 
complexity as it passed from close to the central meridian 
on September 10, to close to the west limb on September 
14. In the photosphere, the total sunspot area fell by close 
to a factor of 2, while the total magnetic flux fell by approx- 
imately a factor of 6. In the corona, the average electron 
density, temperature, and pressure all showed similar de- 
creases in value, which is to be expected considering that 
the plasma in the corona traces out field lines which are ul- 
timately rooted in the photosphere. These results are of con- 
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Figure 10. Top Panel: Plot of Ph versus ^tot for five days of 
MDI and CDS observations. Also shown is a least-squares fit to 
the non-flaring data together with power-laws with slopes of 1 
and 2 from the theoretical predictions. Bottom Panel: Plot of Ph 
versus <I>tot, where Ph has been calculated using non- flaring core 
{dotted line), boundary {dashed line), and quiet-Sun {dot- dashed 
line) parameters from Fig. |6l 



sequence to efforts in understanding active region evolution 
and the relationship between the photosphere and corona of 
solar active regions (e.g., Abbett & Fisher 2003; Ryutova & 
Shine 2004). 

In addition to studying active region evolution, CDS 
and MDI were used to investigate the power-law rela- 
tionships predi c ted b y theoretical models of the corona. 
iMandrini et also investigated theoretical scaling 

laws using magnetic field extrapolations from both vec- 
tor and_Jorigitudinal magnetograms, in light of the work 
of iKlimchuk Porteil (Il995l) . They concluded that mod- 
els involving the gradual stressing of the coronal magnetic 
field are in better agreement with the observational con- 
trains than are wave heating mode ls. This s ame general 
conclusion was also reached by Demo ulin et al. (2003) us- 
ing photospheric and coronal measurements from MDI and 
Yohkoh. Unfortunately these studies relied on broad-band 
filter ratios to determine electron temperatures and den- 
sities; the difficulty associated with making such measure- 
ments is clear from the contradictory results of Priest et al.l 
( 1998). lAschwandenl i200ll) . and lPmlel |200^, who also in- 
vestigated coronal loop heating models using filter ratios 
from Yohkoh/SXT. 

The analysis presented in this paper, on the other hand. 
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is based on well understood line ratio techniques, which of- 
fer a less ambiguous determination of plasma properties. 
With this in mind, the power- law relationship between the 
total heating rate and the total unsigned magnetic flux 
($tot) was detern iined, finding a rel ationship of the form 
^0 81^0.32 l^is^^^ al.^(|l998^ compared the X-ray 
luminosity (which was assumed to be some fraction of the 
total heating power) to active region vector magnetograms, 
finding a similar power- law relation of Lx ^to\^ • The 
result of lFisher et al.l (|l998l) suggests a "universal" relation- 
ship between magnetic flux and the amount of coronal heat- 
ing, regardless of the age or complexity of the active region. 
A similar relationship of Lx ^t^t was foun d using statisti- 
cal samples of late- type stars ( Schriiveiil987l'> . The results of 
this paper therefore lend further observational evidence that 
active regions are heated by magnetically-associated waves, 
rather than multiple nanoflare-type events. 
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